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ABSTRACT 

We present images of infrared emission from the radio jet in 3C66B. Data at three 
wavelengths (4.5, 6.75 and 14.5 /im) were obtained using the Infrared Space Observa- 
tory. The 6.75 /jm image clearly shows an extension aligned with the radio structure. 
The jet was also detected in the 14.5 /xm image, but not at 4.5 /im. The radio-infrared- 
optical spectrum of the jet can be interpreted as synchrotron emission from a popu- 
lation of electrons with a high-energy break of 4xlO n eV. We place upper limits on 
the IR flux from the radio counter-jet. A symmetrical, relativistically beamed twin-jet 
structure is consistent with our results if the jets consist of multiple components. 
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1 INTRODUCTION 

Jets are common in extragalactic radio sources. They are 
believed to trace beams of particles and fields ejected from 
the compact central region at bulk velocities approaching 
the speed of light, combined with material swept up by the 
flow, and are usually visualised through the loss of energy 
as synchrotron radiation. The emission measured at radio 
frequencies should, in the absence of a high-energy cutoff in 
the electron energy spectrum, extend to the optical and in- 
frared bands, although optical counterparts are detected for 
relatively few jets [Martel et al. (1998) list 9 such cases] and 
an even smaller number of infrared jets have been reported. 
One noteable exception is the jet in M87 which has been 
detected in the near-IR (Killeen et al. 1984) and with ISO 
in the mid-IR (Tansley et al. in preparation). 

The radio galaxy 3C 66B has one of the best-studied 
optical jets. The source has a redshift of 0.0212 (Stull et 
al. 1975) and is associated with a 13th magnitude ellipti- 
cal galaxy in a small group in the vicinity of the cluster 
Abell 347. Its radio structure (Leahy, Jagers and Pooley 
1986; Hardcastle et al. 1996) has been studied in detail, 
and reveals a morphology which is intermediate between an 
edge-darkened double and a head-tail structure, possibly in- 
dicating the effect of the parent galaxy's motion through the 
ambient medium. 

The optical jet in 3C 66B was first detected by Butcher, 
van Breugel and Miley (1980). Further ground-based work 
(Fraix-Burnet, Nieto and Poulain 1989) has revealed the 
presence of five distinct knots of polarized optical emission 
along the jet. These correspond to five similar knots in the 
radio emission. A tentative optical detection of the south- 
ern counter-jet was made in the I-band (Fraix-Burnet 1997). 
This is interesting as it allows the radio-optical spectra in 



the jet and the counter-jet to be compared. If the simplest 
twin-relativistic-beam model is valid, then the spectra of 
the jet and counter-jet should differ in intensity and shape 
in ways described by a single component synchrotron model. 

Hubble Space Telescope (HST) images of 3C 66B (Mac- 
chetto et al. 1991) show the optical jet on scales of ~ 0.1 
arcsec, and have revealed an intriguing double-stranded fil- 
amentary structure which resembles the radio image at the 
same resolution. A comparison between HST and VLA data 
at 0.25 arcsec resolution (Jackson et al. 1993) shows further 
similarities. 

Here wepresent data obtained with the Infrared Space 
Observatoryfj (ISO) which appear to detect the jet at 6.75 
and 14.5 ^m. 

Throughout this paper we assume Ho=50 km s" 1 
Mpc -1 . 



2 OBSERVATIONS AND DATA REDUCTION 

3C 66B was observed with the European Space Agency tele- 
scope ISO as part of a larger project to study infrared emis- 
sion from the dusty central regions of a sample of 3CRR 
objects (Tansley et al. in preparation). Observations were 
carried out with the ISOCAM instrument (ISO's near- mid 
IR imaging camera; Kessler et al. 1996, Cesarsky et al. 1996) 
using the LW1, LW2 and LW3 filters to obtain images in 
bands 4.00-5.00, 5.00-8.50 and 12.00-18.00 /im centred on 



ISO is an ESA project with instruments funded by ESA Mem- 
ber States (especially the PI countries: France, Germany, the 
Netherlands and the United Kingdom) and with the participa- 
tion of ISAS and NASA. 
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wavelengths of 4.5, 6.75 and 14.5 /im, respectively. A pixel 
scale of 3 arcsecs was selected, and the angular resolution of 
the data was 3.26, 3.98 and 4.03 arcsecs for the three filters 
(measured from the supplied PSF libraries), LW1, LW2 and 
LW3 respectively. The noise levels on the final images are 
close to the theoretically expected values of 69, 71, 70 ^tJy 
for the three filters, respectively. 

The data were processed initially using the ISOCAM 
Consortium/ESA CAM Interactive Analysis (CIA) package 
(Ott et al. 1997), running under the Interactive Data Lan- 
guage, IDL. The data were corrected for instrumental and 
cosmic-ray glitch artifacts and also for data stabilization us- 
ing the multi-resolution median method.]] The most recent 
(15th March, 1998) calibration-library flat-field and dark 
frames were then applied. The final calibrated images have 
unexpectedly "lumpy" backgrounds, perhaps indicating that 
the supplied flat-fields were not ideal. Further details of the 
analysis are given in Tansley et al. (in preparation) 



3 RESULTS AND ANALYSIS 

The 6.75 /im image can be seen in Figure ^ The contours 
of this image represent the ISO data and the greyscale is a 
1425 MHz VLA image from Hardcastle et al. (1996) 

The final calibrated images were fitted with a three- 
component model, consisting of an extended galactic compo- 
nent (modelled as an elliptical Gaussian based on the optical 
galaxy as a template), an unresolved AGN component, and 
a set of one or two unresolved points representing emission 
from the jet component. The goodness of fit of the model 
at various stages of complexity is shown for the LW2 data 
in Table |l|. Table ^ shows the positions and central offsets 
of the best fit components with this filter. Fitting only an 
extended galactic component or unresolved AGN compo- 
nent gives a significantly worse fit (by A^ 2 =51 and 292 re- 
spectively) than a combination of the two. At this stage, a 
residual map of the ISO LW2 (and LW3) data showed clear 
evidence for excess emission, with structure significantly dif- 
ferent from, and brighter than, the background "lumps". 
This excess was to the east of the AGN/galaxy position in 
the region of the inner radio and optical jet. The overall 
fit improves further (by Ax 2 =14.4 or 22.5) when one or two 
unresolved components representing jet emission were added 
to the model. 

Table |^ presents the fitted flux densities in the three 
ISO filters for the jet components based on the positions 
derived from the LW2 image. We believe the high \ 2 values 
to be caused primarily by inaccuracies in the flat-field cali- 
bration frame rather than the model itself: a \ 2 map of the 
fitted region is similar to a x 2 map of an off-source region 
in showing patchy (scale ~6 - 9 arcsecs) low-level (peaks of 
~0.1 mjy) residual structure. Although this structure lies 
above the nominal noise level, it is much fainter than the jet 
components that we identify: if a post-facto assessment of 
the image noises is made, based on the final % 2 values, the 
errors on the LW1 and LW2 flux densities would increase by 
50 and 68% respectively. Errors in the absolute flux scale can 



t For details see: 


httD: / /isowww. cstcc.esa.nl/manuals/ 


HANDBOOK/III/cam_hb/nodc51.htm: 


: 



also contribute to uncertainties in the measured flux densi- 
ties. The ISOCAM documentation quotes an uncertainty in 
the spectral energy distribution calibration of ±5%. Which 
is much less than the random errors on the jet and is not 
accounted for in our quoted errors. 

Table |j| lists the jet flux densities at radio, IR and op- 
tical wavelengths. The radio measurements were taken from 
images obtained with the VLA at 1.4 and 8 GHz (Hard- 
castle et al. 1996) integrated over the same regions covered 
by the IR components. The corresponding I band flux den- 
sity was taken from the galaxy-subtracted image of Fraix- 
Burnct (1997) obtained at the Canada- France-Hawaii tele- 
scope. The optical flux densities were taken from the HST 
measurements of Jackson et al. (1993), or measured directly 
from archival HST data. However, since the field of view 
(FOV) of the HST Faint Object Camera (FOC) does not 
include the full extent of the IR emission, we assumed that 
the radio-optical spectral index of the entire region was equal 
to that of the observed region and we scaled the measured 
optical flux densities by a factor 1.66 to correct for the 40% 
of the radio emission that lies outside the FOC's FOV. 



4 DISCUSSION 

Figure ^ shows the measurements from Table |i| fitted with 
two separate synchrotron models. In both cases, the volume 
was assumed to be a simple cylinder (of length 5 kpc and 
radius 0.5 kpc) along the length of the IR jet. We assumed 
that the emission arises from a plasma with velocity v=0.7c 
oriented at 45° to the line of sight. These assumptions have 
only a modest effect on the best-fit field, and are consistent 
with the beam decomposition discussed later (Section 4). 

The simplest spectral shape that fits the data is one 
in which p=2 (where N(E)dE = N Q E~ P ) in the radio- 
optical region that turns down sharply at higher energies to 
match the optical points. We fitted two synchrotron spec- 
tra based on indicative models of the population of rela- 
tivistic electrons. The first model population has a high 
energy cutoff at E max =6.5xl0 11 eV and an equipartition 
magnetic field of 3.7 nT. This fit reproduces the optical 
points badly, giving a spectrum that is too steep. The second 
model consists of a broken power law population of electrons, 
where the spectrum breaks from p=2 to p=6 at an energy 
Ei,reafc=4.5x 10 11 eV. The maximum energy in this model is 
E ma:E =15xl0 eV, and the same equipartition field of 3.7 
nT is present. The broken-power-law model shows a consid- 
erably better fit to the optical points (x 2 =215 as compared 
to x 2 — 391 for the model with no power-law break). 

The worst fitting point in both models is the low- 
frequency radio point, though we attribute this to the pres- 
ence of steep-spectrum material in the jet, possibly indica- 
tive of a two-component flow in the jet structure (discussed 
later). Similarly, both fits poorly reproduce the I-band data 
point, though this may merely be a consequence of deficien- 
cies in the simple models since the spectrum near I band is 
strongly dependent on the shape of the break in the electron 
spectrum. Of greater interest is the higher-than-expected 
UV point from the HST observation as this may represent 
an unmodelled component of the electron energy distribu- 
tion. 

To test whether the flatter optical spectrum is an ar- 
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Components Component Flux Density Goodness of Fit 

m Jy X 2 I degrees of freedom 



AGN 


4.4±0.5 


450.1/60 


Galaxy 


8.4±0.4 


208.8/58 


AGN 


2.0±0.3 


157.7/55 


+ Galaxy 


8.2±0.4 




One jet point 


0.4±0.1 


143.3/52 


+ AGN 


1.8±0.3 




+ Galaxy 


8.3±0.4 




Two jet points 


0.8±0.1 


135.2/49 


+ AGN 


0.3±Q.l 




+ Galaxy 


2.1±0.2 






8.2±0.4 





Table 1. Goodness of fit comparison for variations of the model fit at 6.75 (im. 



Component Position of Component Offset From Centre 

RA & Dec of Galaxy/ arcsecs 

Galaxy 02 h 23 m ll?2, 42°59'28" N/A 

AGN 02 h 23 m ll=2, 42°59'28" 

Jet 1 02 h 23 m 12!0, 42°59'35" 13.9 

Jet 2 02 h 23 m ll?6, 42°59'32" 7.2 
Table 2. Positions and offsets of best-fit components. 



tifact of our extrapolation of the observed FOC flux den- 
sities, fluxes were measured on the radio, I-band and HST 
images only within the FOV of the FOC. These data points 
still show a shallower optical spectrum than is predicted by 
our two models. Thus, the extrapolation of the optical flux 
over the full radio-jet region is not the cause of the flatten- 
ing. Furthermore, a U-band observation of 3C 66B using the 
STIS instrument on HST (Sparks et al. in preparation), cov- 
ers the full extent of the radio jet and has a flux density that 
is similar to the scaled Jackson et al. (1993) flux density in 
the same band. 

Forthcoming Chandra X-ray observations of this source 
will allow us to constrain further the model parameters of 
the emission; in particular they will show whether the flat- 
tening of the spectrum in the optical continues to higher 
frequencies. If this is true, it may be the case that we are 
seeing emission in the blue dominated by a newer, flat spec- 
trum component, the origin of which is uncertain. Further 
observations in the optical and UV would be needed to in- 
vestigate this. 

The counter-jet seen in radio maps of Hardcastle et al. 
and possibly in the optical images (Fraix-Burnet 1997) is not 
detected with ISO. Upper limits for the flux densities (Table 
0) have been measured using the radio and I-band maps as 



guides to the emission region. The radio-IR spectral index 
is found to be a > 0.5. 

Fraix-Burnet (1997) identified several components of 
the counter-jet which appeared coincident with radio fea- 
tures. The measured radio and optical flux densities of these 
knots imply that the counter-jet has a flatter radio-to-optical 
spectrum than the jet. This seems to contradict the usual 
interpretation of the jet and counter-jet as physically iden- 
tical, with their brightness asymmetry caused by the rel- 
ativistic Doppler effect. Our IR upper limits lie above the 
fluxes that we expect for the counter-jet, based on a sim- 
ple broken-power-law model with the same parameters as 
the main jet. This model is plotted with the data in Fig- 
ure ^J. As can be seen, the fit to the optical point in the 
counter-jet spectrum is poor, and agrees with the finding of 
Fraix-Burnet (1997) that the spectrum does not conform to 
that expected if 3C66B was a twin beamed structure with 
identical jets. Instead, it is possible that either the source is 
not described by a symmetrical beamed-jet model or that 
the observed optical structures are not real. 

The radio spectrum of the counter-jet spectrum is 
steeper (a = 0.88) than the radio points of the jet (a = 
0.62). A possible explanation could be that we are seeing 
the superposition of an old population of material (forming 
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Filter Goodness of Fit Flux Densities 

X 2 I degrees of freedom Component 1 Component 2 Total 

(mjy) (mjy) (mJy) 

LWlt 110.0/55 <0.36 <0.36 <0.72 

LW2 135.2/49 0.43±0.13 0.67±0.18 1.10±0.22 

LW3 116.0/49 1.08±0.26 0.59±0.18 1.66±0.31 



Table 3. Flux densities for the infrared jet emission, f Jet not detected with this filter. 3<r upper limits are based on a model without 
jet components. Model fitted without jet components. 



Frequency Flux Density Reference 

Hz mjy 

1.30X10 15 0.0061 ± 0.0009 Jackson et al. 1993 (scaled to 

account for smaller FOV) 

9.67X10 14 0.0102 ± 0.0016 Jackson et al. 1993 (scaled to 

account for smaller FOV) 

8.80X10 14 0.0128 ± 0.0014 Macchetto et al. 1991 

7.32 Xl0 14 0.0184 ± 0.0024 This paper (from public archive) 

3.33X10 14 0.111 ± 0.001 Fraix-Burnet 1997 

6.67X10 13 <0.72 This paper 

4.44X10 13 1.10 ± 0.22 This paper 

2.06X10 13 1.66 ± 0.31 This paper 

8.41x10 s 105.0 ± 2.7 Hardcastle et al. 1996 

1.42x10 s 315.0 ± 9.0 Hardcastle et al. 1996 



Table 4. Optical, infrared and radio fluxes for the inner jet. 



a sheath around the jet) combined with newer material re- 
cently ejected and passing through the sheath. Evidence for 
such an enclosing structure has been reported in (for exam- 
ple) Hardcastle et al. (1996), Laing (1996) and Hardcastle et 
al. (1997) and is also described by O'Dea and Owen (1987) 
in relation to the radio maps of NGC 1265. As a weak sheath 
of emission appears at only 1415 MHz in NGC 1265, O'Dea 
and Owen concluded that the sheath has a steeper spectrum 
than the jet itself, as would be expected from older material 
which had undergone spectral ageing. The presence of the 
sheath is attributed to the diffusion of particles out of the 
beam into the static external medium. 

In this model, we expect the inner component to be 
travelling faster than the sheath material, so that emission 
from the inner flow is relativistically-boosted on the jet side 
and dimmed on the counter-jet side, relative to sheath emis- 
sion. Hence the spectrum of the counter-jet would be less 
dominated by the inner material, and consequently steeper. 
Hardcastle et al. (1996) have produced a polarization inten- 
sity map that shows clearly a distinct line of depolarization 



marking the boundary between the sheath and the inner 
jet in 3C 66B. Flux-density measurements of the inner and 
outer regions of the jet show that the inner jet is dimmed 
on the counter-jet side by a greater amount than the sheath 
component, as the model predicts. 

If we accept this model of the inner beam and outer 
sheath then we can use the measured flux densities from the 
radio polarization intensity map to estimate the inner flow 
and sheath velocities. We can represent the emission from 
the jet as the sum of emission components. 

q T/ r a s heath+2 . rfl inn er + 2 

Stotal — ^-sheath v sheathO sheath dinner V inner ^ nner 

where <5 is the bulk relativistic Doppler factor: 

5 = j' 1 (1 - pcosey 1 

e is the emissivity of the material, V is the volume of the 
emitting region and a is the spectral index of the region 
(either sheath or inner). 6 is the orientation angle to the 
line of sight and /3 is the ratio of the velocity to the speed of 
light, c. Measuring equal volumes on both sides and equating 
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5 10 15 20 5 10 15 20 




02 23 14 13 12 11 10 02 23 14 13 12 11 10 

RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) 



Figure 1. Image showing the ISO LW2 data, with galaxy emission included (left) and subtracted (right), as contours (lowest contour 
at 3o-; levels O.llmjy/pixel X -1, 1, 2, 3, 4, 5), with VLA image at 1425 MHz (Hardcastle et al. 1996) shown in grey-scale (0 to 20 
mjy/beam). Note the disappearence of the companion galaxy from the subtracted image (right) 



5 10 15 20 5 10 15 20 




02 23 14 13 12 11 10 09 02 2314 13 12 11 10 09 

RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) 



Figure 2. Image showing the ISO LW1 (left) and LW3 (right) data as contours (lowest contour at 3<r [0.09mjy/pixel and 0.18mjy/pixel 
X -1, 1, 2, 3, 4, 5 respectively] level), with VLA image at 1425 MHz (Hardcastle et al. 1996) shown in grey-scale (0 to 20 mjy/beam). 
The companion galaxy can be seen clearly at 02 h 23 m 13!0, 42°59'17'f 



each component, we find Pinner ~ 0.7 and /3 s heath ~ 0.2 
(using ai nner — 0.62 and a s h ea th = 0.88 measured from the 
jet and counter-jet radio spectral indices respectively and an 
orientation angle to the line of sight, 8 = 45°). 

The observed steepness of the spectrum of the counter- 
jet in the radio initially places doubt on the association of 
the I-band detections with material in the jet itself. How- 
ever, if the multi-component model proposed here is ac- 



cepted, then at higher frequencies the counter-jet may be 
dominated by the inner jet rather than the steeper spec- 
trum sheath, despite the relativistic dimming of the inner 
jet. Using this assumption, we can perform a simple cal- 
culation to test whether the optical flux of the counter-jet 
from the Fraix-Burnet image is what we expect to observe. 
In fact, using the Doppler parameters of the inner jet, the 
optical spectral index of a—lA (Jackson et al. 1993) and 
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Frequency (Hz) Frequency (Hz) 



Figure 3. Optical-IR-Radio spectrum of 3C 66B showing fitted synchrotron profile. Top: calculated using an electron spectrum with a 
high-energy cutoff. Bottom: calculated using an electron spectrum with a high-energy break. In both cases the frequencies and the flux 
densities are in the rest frame of the inner-jet component, taking /3=0.7 and 6>=45° (see Section 4). 



the flux density measurement from the corresponding fre- 
quency in the jet, we find an expected flux density of 2 /ijy, 
much lower than the flux of 7 /xJy measured by Fraix-Burnet 
(1997). 



5 CONCLUSIONS 

We have made the first mid-IR detection of emission from a 
radio jet. Fits to synchrotron models using radio and opti- 
cal data points have shown that the overall spectral-energy 
distribution of 3C 66B resembles a simple synchrotron spec- 
trum from the radio to the optical, though the optical slope 
is unexpectedly flat, perhaps because of inhomogeneities in 
the jet. The counter-jet in 3C66B is not detected in our 
data. The jet and counter-jet spectra can be interpreted in 
terms of a model in which the beam structure is composed 
of a slow-moving (/3=0.2), steep-spectrum component and a 
faster-moving ((3=0.7), flatter-spectrum component, a con- 
clusion supported by radio images of polarised emission from 
the jet structure. 



Frequency 
Hz 


Flux Density 
mjy 


References 




3.33X10 14 


0.007 


Measured from the maps 


of 






Fraix-Burnet 1997 




6.67X10 13 


<0.16 


This paper. 




4.44 xlO 13 


<0.15 


This paper. 




2.06X10 13 


<0.29 


This paper. 




8.41 xlO 9 


15 


Measured from the maps 


of 






Hardcastle et al. 1996 




1.42X10 9 


73 


Measured from the maps 


of 






Hardcastle et al. 1996 





Table 5. Counter-jet infrared and radio fluxes. 



by the ISOCAM PI, C. Cesarsky, Direction des Sciences de 
la Matiere, C.E.A., France. 
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Figure 4. Radio-IR-Optical spectrum of the counter-jet in 3C66B showing fitted high-energy cutoff synchrotron profile. Rest-frame 
frequencies shown. The frequencies and flux densities shown are in the rest frame of the inner-jet component, taking /3=0.7 and 6=45° 
(see Section 4). 
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